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NONSINGULAR EFFICIENT MODELING OF 
ROTATIONS IN 3-SPACE USING THREE COMPONENTS 

NORMAN J. GOLDSTEIN 

o 

CN ' Abstract. This article introduces yet another representation of rotations in 

►^,^' 3-space. The rotations form a 3-diinensional projective space, which fact has 

rrt ■ not been exploited in Computer Science. We use the four afhne patches of this 

projective space to parametrize the rotations. This afhne patch representation 
is more compact than quaternions (which require 4 components for calcula- 
tions), encompasses the entire rotation group without singularities (unlike the 
Euler angles and rotation vector approaches) , and requires only ratios of linear 
or quadratic polynomials for basic computations (unlike the Euler angles and 
rotation vector approaches which require transcendental functions). 
Tf^ ' As an example, we derive the differential equation for the integration of 

angular velocity using this affine patch representation of rotations. We re- 
mark that the complexity of this equation is the same as the corresponding 
1 ' quaternion equation, but has advantages over the quaternion approach e.g. 

f_\ ' renormalization to unit length is not required, and state space has no "dead" 

directions. 



^ ■ 1. Introduction 

^ . . 

^^ ' This work was motivated by the need to have a practical and robust way to 

ly-v , model rotations and projective spaces in the software that we write. These are 

f^ ' applications in photogrammetry, remote sensing and simultaneous localization and 

^^ , mapping. Rotations are used in modeling the poses of the cameras, while the 

focal planes of the cameras are projective spaces. Although it is true that the 
imaging focal plane is contained in a single afHne patch of RP"^, relevant vanishing 
L^ ■ points may lie on, or near, the line at infinity, so that the coordinates of the imaging 

^ ^N , affine patch are not appropriate for robustly describing and managing the vanishing 

C^ ' points. 

Traditionally, (unit) quaternions have been used to model rotations, as the 
quaternions have only 4 components - just one more parameter than the dimension 
of the space of rotations (see, for example [U [2]). However, the extra component 
is bothersome, requiring extra storage, requiring extra work to manage the redun- 
dancy (renormalization to unit length), and likely affecting the convergence rate 
due to motion in the redundant length direction of the quaternion, when we are 
optimizing to find the best rotation for a particular problem. 

Another approach to representing rotations is to use the three Eulerian angles^ 
e.g. [3]. However, this approach suffers from not being everywhere invertible, so 
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is not a global solution to representing the rotations. Moreover, the transcenden- 
tal transformation between Eulerian angles and rotation matrices is not ideal for 
software implementation. A similar situation is to be found with the rotation vec- 
tor representation of a rotation. The on-line article [4j discusses this and other 
representations of rotations. 

In the following sections, we review the standard coordinates on projective space, 
and then derive the differential equation for the integration of angular velocity using 
the projective space manifold structure (afhne patches) on the rotation group. We 
will see that with this novel representation of the rotations, the differential equation 
has the same complexity as the corresponding quaternion differential equation, and 
has other desirable properties over the quaternion case. 

2. Notation 

All vector spaces and projective spaces are over the real numbers. For a point 
z in n-dimensional projective space, P", its homogeneous coordinates are denoted 
z= [z° , z^ , . . . z"-]. 

The traditional manifold structure on P" is given as 

• The k*^^ affine patch, fc = 0, . . . n, on P" is defined to be the set 

Vk^{z^[z°,z\...z"]:z''^0} 

• The coordinate chart for Vk is 

defined by the formula Xk{z) = zj z^ and throw away the 1 in the /c*'* slot. 
This is seen to be a 1-1 map onto all of M". 

The following algorithm template illustrates how to use the manifold structure 
on P" to iteratively step along a sequence of points in projective n-space. For 
simplicity, we do not specify what it means to "iterate" , and we omit any stopping 
criteria. 

Algorithm 2.1. Iterate on Projective Space 

Input An initial point, zo £ P"- 

(1) Let k be the index of the largest component, |zo'^|- 

(2) Represent zq in the fc*'' affine patch. So, initially, all the coordinate com- 
ponents in this patch have magnitude at most 1. 

(3) Iterate to the next value of z, in this patch. 

(4) If all coordinate components of z have magnitude at most 2, go back to 
Step 3. Otherwise, proceed to Step 5. 

(5) Reconstruct the homogeneous coordinates, z, for the current point, and let 
k be the index of the the largest component \z^\- 

(6) Represent z in the A;*'' affine patch, and proceed to Step 3. 

In Step 6, some book-keeping may be required to keep the overall environment of 
the algorithm consistent when switching patches. The choice of the value, 2, in 
Step 4 is somewhat arbitrary, but has worked well for us. The point is to monitor 
when a particular patch is no longer appropriate for the current value of z, and 
then to switch to a good patch. If z begins to converge, the switching of patches 
will cease. 
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This algorithm may be used when optimizing for a "best" value of z, or when 
incrementally solving a differential equation on projective space. In the next section, 
we write down explicitly the differential equation for integrating angular velocity, 
using the projective space manifold structure of the rotations. 

3. Integrating Angular Velocity 

We use the convention that the components of a quaternion, q, are indexed as 
q^,q^,q^ and q^, with s = q'^ being the scalar component, and v denoting the 
3-vector of the other (spatial) components. 

The particular rotation matrix associated to q is 

R = /(2s2 - 1) + 2w* + 2sMy 
Where, / is the identity 3x3 matrix, and Ma is the usual 3x3 skew matrix defined 

by 

Mab — a X b 
Explicitly, 

Ma = 

Our jumping off point for analyzing angular velocity is the well-known ([l][2]) 
differential equation relating angular velocity, uj, and quaternions: 

This corresponds to the matrix equation 

at 
It is Equation 13. II that we wish to translate to the four affine patches of the rota- 
tions. To that end, we make explicit the components of the matrix in Equation l3.1l 
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(3.2) H = 



/ -c^i -^2 _^3 \ 

Cjl w3 -^2 

a;2 -cj^ cji 



We restrict attention to patch i, i = 0,1,2,3, and let 



X 
^2 



be the 3 components parametrizing the rotations from this patch. The correspond- 
ing homogeneous coordinates, x are obtained by forming the 4- vector with a 1 in 
the i*^ slot. So, a; is a non-unit quaternion, and we define a time-varying scale 
factor, a = a(t) by the equation 

X = aq 
where q is now a unit quaternion. Consequently, Equation 13.11 is valid for q, and 
we may write 

dx da 1 

— = — q + a—Hq 
dt dt^ 2 ^ 
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which may be rewritten as 

(3.3) t'^'^l"' 

where f3 = -^ I ol. To reduce this equation, further, we define the 3x4 matrix whose 
columns are the non-zero entries from the corresponding columns of H. 



^1 


-^1 


-u? 


-^3 


u? 


-U^ 


W3 


-W2 


L.3 


W2 


-L.1 


c^i 



(3.4) W = 

We recall that the i*'' component of x is identically 1, so its derivative is 0, and we 
conclude from the i*'' component of Equation 13.31 that 

where iJ* is the z*'' row of H . This can be rewritten as 

(3.5) ^=\w,-x 

where Wi is the i*'' column of W. 

Extracting, now, the non-i*'' components of Equation 13.31 we obtain 

(3.6) ^ = /3a; + i(Ty,+ijWo 

where H^^ is the 3x3 matrix obtained from H by deleting the z*'' row and the i*'' 
column. By inspecting the four cases, we see that 

(3.7) il« = {-Xf^^Mw^ 
Inserting Equations 13.51 and 13.71 into Equation 13.61 we obtain 

dx \ 1 

^ - ^{W,-x)x^^{W,^{-\y+^W,y.x) 

(3.8) = i[VK, + (iy, •a;):r + (-l)'+iW^, xx] 



Equation 13.81 is the 3-dimensional differential equation to integrate angular ve- 
locity using affine patches on the rotation group. The RHS contains 12 multiplica- 
tions, excluding the scaling by 1/2. This is the same number of multiplications as 
in Equation 13.11 which is the corresponding quaternion differential equation. 



3.1. Difference Equation. In practice, a difference equation is used to integrate 
angular velocity. This section shows how to obtain identical results to the quater- 
nion solution, but by using the affine patch representation of the rotations. This 
has the advantage that explicit renormalization to unit length is not required, as it 
is inherent to the affine patch approach. 

For a time step, of length Ai, we wish to find the x increment. Ax, so that the 
two expressions, q -f ^Hq/^t and x + Aa:, represent the same rotation i.e. 

{q + \HqMy 

where i denotes all the components except for the i*''. Divide the top and bottom 
of the LHS by g* to obtain 
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(3-9) ^^TTTTTT = ^ + ^^^ 



l + ^H^xAt At' 



We recall from the derivation of Equation 13.81 that 

Hh = W,- (-l)'Wi X X 

Wx = -W^ ■ X 

so that Equation 13.91 simplifies to 

. . Ax _ W^ + i^iy+'W^ xx + {W,- x)x 

^ ' ^ At 2-{W^■ x)At 

It is interesting to note that this value of Ax/ At is the same as dividing dx/dt 
by the scale factor 1 — {Wi ■ x)At/2. 

4. Conclusions 

We have applied elementary manifold theory to construct a global, efficient, 
singularity-free parametrization for the rotation group. We have derived differential 
and difference equations for integrating angular velocity using this novel representa- 
tion of the rotations, and have shown that the complexity of using these equations 
is on a par with the corresponding quaternion equations. On a theoretical basis, 
however, the difference equation derived in this article is more sound, as there are 
no renormalizations to unit length required. 
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